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ABSTRACT 

We present results of our Power Spectral Density (PSD) analysis of 30 AGN using the 58 month 
light curves from Swift's Burst Alert Telescope (BAT) in the 14-150 keV band. PSDs were fit using 
a Monte Carlo based algorithm to take into account windowing effects and measurement error. All 
but one source were found to be fit very well using an unbroken power law with a slope of ~ —1, 
consistent at low frequencies with previous studies in the 2-10 keV band, with no evidence of a break 
in the PSD. For 5 of the highest S/N sources we tested the energy dependence of the PSD and found 
no significant difference in the PSD at different energies. Unlike previous studies of X-ray variability 
in AGN, we do not find any significant correlations between the hard X-ray variability and different 
properties of the AGN including luminosity and black hole mass. The lack of break frequencies and 
correlations seem to indicate that AGN are similar to the high state of Galactic Black Holes. 
Subject headings: galaxies: active - galaxies: nuclei - galaxies: Seyfert - X-rays: galaxies 



1. INTRODUCTION 

Active Galactic Nuclei (AGN) are thought to be 
powered by accretion onto a Supermassive Black Hole 
(SMBH) because of their high luminosity and small size. 
The gravitational potential energy of infalling material 
in the accretion disk is converted into thermal energy 
and radiated away. However infall only explains AGN 
spectra up to UV energies whereas AGN are also seen 
to emit strongly in X-rays. To explain this, a corona 
must be invoked that is thought to Compton upscatter 
the UV photons from the inner accretion disk to X-ray 
energies. The upscattered photons are radiated quasi- 
isotropically so some directly travel to the observer or il- 
luminate the accretion disk. This illumination produces 
the so-called "Compton hump" seen in the spectra of 
AGN that peaks around 30 keV (eg. George & Fabian 
(1991); Nandra & Pounds (1994) for a theoretical dis- 
cussion and Magdziarz & Zdziarski (1995); Rivers et al. 
(2011) for observational evidence.) 

AGN X-ray emission also exhibits very strong and 
rapid variability that suggests the origin of the emission 
is very near the SMBH itself. By studying the spec- 
tral and temporal properties of the X-rays, we can learn 
about the environment inside the strongest gravitational 
potential wells in the universe. 

Variability has been a known property of AGN for sev- 
eral decades. Barr & Mushotzky (1986) found that for 
Seyfert galaxies the timescalc of X-ray variability is in- 
versely related to the X-ray luminosity. Using EXOSAT, 
Lawrence et al. (1987); Lawrence & Papadakis (1993); 
Green et al. (1993) were the first to be able to charac- 
terize the variability as aperiodic following a power law 
relationship with frequency, -P(i^) oc i^~", with P(i^) be- 
ing the power at a given frequency and a ~ 2. 

The most powerful and often used tool in study- 
ing time variability is the Power Spectral Density 
function (PSD), which measures the contribution to 
the total variance that each temporal frequency pro- 
vides. Papadakis & McHardy (1995) created the first 
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broadband PSD with EXOSAT observations while 
Edelson & Nandra (1999), using the new i?ossi X-ray 
Timing Explorer (RXTE) that provided well sampled X- 
ray light curves, constructed the first high quality broad- 
band PSD and measured the first "break" in the PSD 
where a changed from a value of 1.74 to 0.73 as tem- 
poral frequency decreases. The slopes measured were 
later found to be consistent with 2 and 1 respectively due 
to Edelson & Nandra (1999) not accounting for aliasing 
and a miscalculation of errors on the PSD (Uttley et al. 
2002). This break was theoretically needed or else the 
variability would diverge as timcscales increased. 

Since Edelson & Nandra (1999) there have been 
numerous PSD studies of AGN (Uttley et al. 2002; 
Uttley & McHardy 2005; McHardy et al. 2004, 2005; 
Markowitz et al. 2003b) all aimed at studying the long 
and short term variability of AGN, measuring the break 
frequencies and power law slopes of the PSD. From 
all these studies, a relationship emerged between the 
break frequency and black hole mass, however not only 
for AGN but including Galactic Black Holes (GBH) 
(Papadakis 2004; Done & Gierlihski 2005) . Since AGN 
and GBH both harbor black holes at their centers, it 
had long been postulated that AGN are just scaled up 
versions of the same system with the same processes gov- 
erning both. Break frequencies also occur in the PSDs 
of GBHs and using simple scaling with black hole mass, 
a rough correlation was found leading to speculation the 
break frequency is associated with the inner edge of the 
accretion disk. However there was a large scatter in the 
correlation indicating that a second variable was needed 
to strengthen relationship. In particular it was seen that 
at the same black hole mass, higher accretion rate Nar- 
row Line Seyfert I's (NLSl) displayed a higher break fre- 
quency than "normal" Seyfert Is, so accretion rate (ac- 
tually bolomctric luminosity to make all parameters in- 
dependent) was introduced into the relationship. These 
two properties of black hole mass and accretion rate cor- 
rectly account for the scaling of break frequency from 
stellar mass X-ray binaries up to supermassive AGN. 
(McHardy et al. 2006). 
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All of these PSD studies, however focus on the 2-10 
keV energy range because that was the best band for 
RXTE. This band though can be affected by absorp- 
tion, so most of the samples have been limited to Seyfert 
I's since these, according to the unified model, have a 
clear view of the accretion disk. Harder energies with 
E > 10 keV, though, "see through" any gas as long 
as it is "Compton thin" with hydrogen column density, 
A'^H < 3 X 10^"* cm"^ With the Swift satellite and its 
Burst Alert Telescope (BAT) , we can now begin to study 
the hard X-ray sky in an unbiased way because its energy 
range is 14-195 keV. BAT has been continuously moni- 
toring the sky with a wide field of view for ^6 years, thus 
providing high quality light curves for many sources for 
the use of PSD analysis. 

Since Compton refiection peaks in this energy 
range we can also begin to answer the question of 
whether the continuum and reflection components of 
AGN spectra show different variability characteristics. 
Sobolewska & Papadakis (2009) analyzed RXTE spectra 
of 10 AGN and found that the model that best described 
the variability was one with a constant refiection compo- 
nent and an intrinsic power law that varies in both flux 
and power law index. Miniutti et al. (2007) studied the 
Seyfert 1 galaxy, MCG-6-30-15 and also found reduced 
variability in the 14-45 keV band where reflection peaks 
in the energy spectrum. 

Beckmann et al. (2007) looked at the hard X-ray vari- 
ability using the 9 month BAT catalog, however the 9 
month light curves weren't of high enough quality for 
a PSD study. Instead two other tools, excess variance 
and structure function were used in their analysis. They 
found that absorbed, Seyfert 2 galaxies show more vari- 
ability than unabsorbed, Seyfert 1 galaxies, which is sur- 
prising since unless the AGN is Compton thick there 
should not be a large effect on the variability espe- 
cially at these long timescales. They also reported an 
anticorrelation between luminosity and variability, sim- 
ilar to that seen in softer X-rays suggesting the same 
variability process is working in the 2-10 keV energies 
as in 14-195 keV energies. However recent work by 
Emmanoulopoulos et al. (2010) have highlighted the un- 
reliability of structure functions in determining variabil- 
ity properties. 

With the 58 month light curves now available, we con- 
struct and measure the PSDs of a sample of BAT se- 
lected AGN to further study the hard X-ray variability. 
We provide the very first PSD slopes and normalizations 
in the 14-150 keV energy band and compare the PSD 
measured variability with AGN properties to determine 
correlations. Because of the large energy range of BAT, 
we also are able to test any underlying energy dependence 
of the PSD as suggested by Markowitz et al. (2007) for 
Mrk 766 and Cui et al. (1997)for Cygnus X-1. 

The paper is organized as follows. §2 details the BAT 
observations and reduction of light curves including re- 
binning and filtering. §3 describes the construction of 
PSDs, selection of sample, and determination of PSD 
parameters. §4 highlights the results of the model fitting 
including the energy dependence of the PSD and correla- 
tions between AGN properties and variance. §5 discusses 
the implications of these results and compares them to 
previous studies of AGN variability. 



2. DATA REDUCTION 

Since November 2004, Swtft (Gchrcls et al. 2004) with 
the Burst Alert Telescope (BAT) (Barthelmy et al. 2005) 
has been searching for gamma-ray bursts (GRBs) in the 
hard X-ray sky (14-^195 keV). Swift/BAT uses a 5200 cm^ 
coded-aperture mask above an array of 32768 CdZnTe 
detectors to produce a wide field of view ( ^ 1 str fully 
coded, 3 sr partially coded) of the sky. When a GRB 
appears. BAT software can quickly locate the position 
to within 1-4 arcmin using a transformation algorithm 
and allow other instruments and telescopes to observe 
the GRB. 

However while in survey mode (ie not specifically ob- 
serving a GRB), BAT is continuously observing the sky. 
Every five minutes a new snapshot image is produced, 
and these series of snapshot images are then cleaned so 
count rates can be extracted for sources that were within 
the field of view. In this way, snapshot images are turned 
into light curves for a large number of hard X-ray sources 
due to the wide field of view and large sky coverage of 
BAT. (See TucUcr et al. (2010) for details of cleaning 
process.) Each lightcurve contains count rates over a 
^5 year period in eight different energy bands: 14-20, 
20-24, 24-35, 35-50, 50-75, 75-100, 100-150, 150-195 
keV, and then a total count rate for the entire 14-195 
keV band with each datapoint corresponding to single 5 
minute snapshot image. 

For our PSD analysis, we first downloaded the 58- 
month raw snapshot lightcurves for each AGN provided 
by Baumgartner (2012) on HEASARC^ We then added 
together the count rates of the first 7 bands, producing 
a 14-150 keV band light curve to maximize S/N in the 
PSD. We decided to discard the 150-195 keV band data, 
because we felt it did not add any information to the 
variability of the sources due to the low S/N of the light 
curve at these high energies.^ The raw lightcurves were 
rebinned into time bins with a width of 5 days using a 
weighted average so that the PSD was above white noise 
level that overwhelms at high temporal frequencies due 
to Poisson noise. With this rebinning, our PSDs span 
the frequency range 10~^ — 10~^'^ Hz, 2 orders of mag- 
nitude. Also any bin in the light curve with a fractional 
exposure < 2% was not used in the analysis to ensure 
Gaussian statistics were still a good assumption. A cut 
of < 2% was found to routinely keep the number of raw 
light curve points in each bin above 20, but still allowed 
for some bins to be kept that had less than 20 but a 
large exposure time. Gaps in the light curves due to the 
filtering process. South Atlantic Anomaly, GRB obser- 
vations, and satellite maintenance, were filled by linearly 
interpolating between the edges of the gaps without sig- 
nificantly altering the PSD. If too many bins are missing 
in the light curve, linear interpolation could cause an 
artificial smoothing and reduce the variability at high 
frequencies. All of the sources in our analysis had less 
than ^20% of their light curve missing due to gaps. Fig- 
ure 1 shows the light curves of a high S/N source (Cen- 
taurus A) and a source 10 times fainter (NGC 3516) af- 
ter all rebinning, filtering, and interpolation have taken 



^ http : //heasarc . nasa . gov/docs/swif t/results/bs58mon/ 
^ Private communication with Jack Tucllcr and Wayne Baum- 
gartner 
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place. Notice there is still distinct variability even at long 
timescales. 

3. DATA ANALYSIS 

3.1. Power Spectral Density Function 

To construct the PSD, we first subtracted the mean 
count rate, /x, from the light curve, which has the effect 
of removing the zero frequency power from the PSD. We 
then estimate the PSD using the periodogram which is 
the squared modulus of the Discrete Fourier Transform 
(Oppenheim & Schafer 1975), 
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where Vj = y , j = 1 . . . ^y , T is total duration of the 
light curve, N is the number of bins in the light curve, 
and Xi is the count rate at time ti. The periodogram is 
then defined as 
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With this normalization, —2772 , known as fractional rms- 
squared normalization, periodograms from different seg- 
ments of the light curve can be equally compared as 
well as PSDs from different AGN. Also the integral 
of the periodogram over a frequency range is equal 
to the contribution to the fractional rms-squared vari- 
ance due to variations from the corresponding timescales 
(Miyamoto et al. 1991; van der Klis 1997). However 
variability is also introduced due to measurement errors 
and will manifest itself as white noise (equal amplitude 
at all frequencies) in the PSD and cause flattening of the 
observed PSD at high frequencies. The expected noise 
level in the PSD from a light curve with Gaussian errors 
is given by (Vaughan et al. 2003a). 
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'^orr — Si "^err i where (Terr.i IS the crror on the count rate 
at time ti. To test this estimate of the noise level, we 
constructed PSDs of a few clusters of galaxies. Clusters 
are constant sources of hard X-ray emission (Wik et al. 
2011) so the PSD should be fairly flat at the Pnoise level. 
Figure 2 shows the PSD of the Coma cluster, one of the 
brightest galaxy clusters in the BAT catalog. As can be 
seen, Coma's PSD (solid line with error bars) is approx- 
imately flat at the noise level (dashed line) within the 
error bars that were calculated using the method from 
Papadakis & Lawrence (1993). Therefore Equation 3 
should be a good estimate of the white noise variabil- 
ity power from measurement errors. 

Following Papadakis & Lawrence (1993), we then loga- 
rithmically binned the periodogram every factor of 1.6 in 
frequency creating equally spaced bins of width ^0.2 in 
logarithm. The lowest two frequency bins however were 
widened so there were at least 3 original raw PSD points 
used in the average. As Papadakis & Lawrence (1993) 
also show, a bias of ^-.253 is introduced into the power 
of the periodogram so this is corrected for at this stage. 



To determine which sources to use in our model fitting, 
we constructed the observed PSDs of all the AGN in the 
BAT catalog that had a S/N greater than 20 (76 sources) 
using Equation 2 and used visual inspection to look for 
a significant detection of power in the PSD well above 
Pnoise- S/N was taken from the BAT catalog and is the 
significance of source detection in the BAT sky map. 30 
sources showed a significant detection of power and thus 
qualified for further analysis with PSD fitting and are 
listed in Table 1 along with the properties of the BAT 
observation over the length of the mission. All signal - 
to-noise ratios and AGN classifications were taken from 
the BAT 58 month catalog'^. Mean count rates were de- 
termined by arithmetically averaging all the count rates 
together after rebinning and filtering. Missing fraction 
was determined by taking the ratio of the number of 
time bins in the rebinned light curve to the number of 
bins expected if the light curve was uniformly and fully 
sampled. 

In total our sample consists of 3 Seyfert I's, 11 Seyfert 
2's, 9 Seyfert 1.5's, 1 Seyfert 1.9, 3 Blazars, 2 NLRG's, 
and 1 BLRG. Because our sample was selected through 
visual inspection, there is an inherent bias towards bright 
and/or highly variable objects. Blazars are known to 
be highly variable, much more than Seyfcrts so these 
were expected. But what is interesting is many more 
Seyfert 2's than Seyfert I's were selected. Our original 
76 source sample from the catalog contained 33% Seyfert 
2's (25/76), 28% Seyfert I's (21/76), and 20% Seyfert 
1.5's whereas the sample we are fitting has 36% Seyfert 
2's, 10% Seyfert I's, and 30% Seyfert 1.5's. However it 
seems that the sample we have selected is mainly those 
sources that have a high S/N so the distribution differ- 
ence is more of an indication that there are more high 
S/N Seyfert 2's than there are Seyfert I's. Figure 3 shows 
the distribution of S/N for our selected sample and those 
that were not selected. We clearly have selected those 
objects that have a high S/N because these objects are 
the ones that will most likely have enough intrinsic vari- 
ability above that of noise assuming that there is not a 
radical difference in intrinsic source variance. Figure 3 
also shows the luminosity distribution between the se- 
lected and unselected sample. It seems we are not biased 
towards high or low luminosity which means if there ex- 
ists a correlation between the variance and luminosity we 
should be able to observe it. 

Table 2 contains properties of each AGN determined 
from previous studies. X-ray luminosity was taken from 
HE AS ARC. Bolometric luminosity was calculated using 
the strong correlation between Lboi and Lx empirically 
found by Winter et al. (2012). 



log(Lboi) = 1.1157 log(Lx) - 4.228 



(4) 



This relation however was only found for Seyfert 1-1.5 so 
it is possible that for Seyfert 2's this is an underestimate 
of the bolometric luminosity since they have a higher hy- 
drogen column density which absorbs the X-ray emission. 
However for the BAT band, iVn needs to be > 10^^ cm^ 
to have a significant effect. In our sample only 2 sources, 
NGC 4945 and Mrk 3, have a column density this high so 
we don't expect a serious discrepancy between the true 
bolometric luminosity and our estimate. 

http : //heasarc . nasa . gov/docs/swif t/results/bs58mon/ 
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Figure 1. Swift/HAT light curves of Ccntaurus A and NGC 3516 after rebinning, filtering, interpolation. Interpolated points arc shown 
in red without error bars. 

used the standard relation Lodd = 1-26 x IO'^^Mbh- 

Figure 4 shows the observed PSDs for two of our 
brightest AGN, Cen A and NGC 4151 as well as their 
estimated Poisson level from equation 3. Most of the ob- 
served PSDs have a characteristie downward slope, but 
there is varying normalization. To quantify the best-fit 
parameters we use a teehnique based on Monte Carlo 
simulations that is described in the next section. 

3.2. Model Fitting 

There are three main complications that arise when 
attempting to fit a model to an observed PSD. The first 
two are windowing effects known as red-noise leak and 
aliasing. 

Red-noise leak is the transfer of power from frequencies 
lower than the minimum frequency sampled by the ob- 
served light curve (i.e. T^^, Tobs is the length of the light 
curve) to higher frequencies. This occurs due to the finite 
length of the light curves, where we have essentially mul- 
tiplied the intrinsic, infinitely long light curve by a box 
function with a width equal to the length of the observed 
light curve. 

Aliasing is the transfer of power from frequencies 
higher than the maximum frequency sampled by the re- 




log(,.) (Hz) 

Figure 2. Observed PSD of the Coma galaxy cluster (solid line 
with error bars) plotted with the estimation of the Poisson noise 
contribution using Equation 3. 



Black hole masses were all taken from the literature as 
well as A^H, the intrinsic neutral hydrogen column den- 
sity. The references for both Mbh and A^h are listed in 
the table along with the method used to determine each 
black hole mass. For Ledd, the Eddington luminosity, we 
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Table 1 

AGN Observation Properties 









Mean Count Rate 


Fraction 


Name 


Type 


S/N 


(cts s-^) 


Missing 


(1) 


(2) 


(3) 


(4) 


(5) 


30 273 


Blazar 


156.8 


6.38E-04 


0.17 


30 390.3 


Syl 


47.97 


1.64E-04 


0.07 


30 454.3 


Blazar 


40.43 


1.83E-04 


0.15 


4U 1344-60 


Sy 1.5 


38.73 


1.66E-04 


0.12 


Oen A 


NLRG 


428.7 


2.08E-03 


0.16 


Oircinus Oalaxy 


Sy 2 


101.7 


4.62E-04 


0.10 


Oygnus A 


NLRG 


54 


2.14E-04 


0.08 


ESO 506-G027 


Sy 2 


26.44 


1.38E-04 


0.17 


10 4329A 


Syl 


101.1 


5.03E-04 


0.18 


IGR J21247+5058 


BLRG 


83.3 


3.27E-04 


0.08 


MOG-05-23-016 


Sy 1.9 


90.4 


3.55E-04 


0.13 


MOG+08-11-011 


Sy 1.5 


49 


2.31E-04 


0.19 


Mrk 110 


Syl 


30.42 


9.49E-05 


0.10 


Mrk3 


Sy 2 


48.5 


2.12E-04 


0.07 


Mrk 348 


Sy2 


70.4 


2.58E-04 


0.13 


Mrk 421 


Blazar 


109.5 


3.10E-04 


0.10 


Mrk 6 


Sy 1.5 


26.75 


9.02E-05 


0.07 


Mrk 926 


Sy 1.5 


45.57 


1.80E-04 


0.19 


NOG 2110 


Sy 2 


98.1 


4.98E-04 


0.13 


NOG 3227 


Sy 1.5 


56.2 


1.73E-04 


0.14 


NOG 3516 


Sy 1.5 


62.3 


1.76E-04 


0.06 


NOG 3783 


Sy 1.5 


68.7 


2.66E-04 


0.16 


NGG 4151 


Sy 1.5 


275 


8.37E-04 


0.07 


NOG 4388 


Sy2 


110.7 


4.22E-04 


0.12 


NOG 4945 


Sy 2 


76.1 


3.72E-04 


0.13 


NGG 5252 


Sy2 


42.4 


1.81E-04 


0.15 


NGG 5548 


Sy 1.5 


32.55 


1.27E-04 


0.07 


NOG 6814 


Sy 2 


22.52 


1.14E-04 


0.18 


NGG 4507 


Sy2 


64.6 


2.82E-04 


0.16 


NOG 7172 


Sy 2 


60.1 


2.48E-04 


0.21 



Note. — Col.(l): Name of AGN. Col. (2): Type of AGN. Col. (3): 
Signal-to-noise ratio. Col. (4): Mean count rate over the whole light curve. 
Col. (5): Fraction of missing data due to gaps in the light curve. 



binned observed light curve (i.e. (2Ari-obin)~"'^ or the 
Nyquist frequency, AT,cbin is the width of each time bin 
of the hght curve) to lower frequencies and occurs due to 
the discrete sampling of the observed light curve. This 
can be thought of again as multiplying the continuously 
sainpled intrinsic light curve by a series of box functions 
that have a width equal to the sampling time. Both of 
these effects multiply the true light curve by a window 
function in the time domain which is manifested as a 
convolution in the Fourier domain. Both are also model 
dependent effects, because the strength of each is depen- 
dent on the true variability power in the corresponding 
frequencies for each effect. A priori we have no prescrip- 
tion for correcting for them unless we assume a PSD 
shape, (see Papadakis & Lawrence 1993; Vaughan et al. 
2003a, for a more detailed discussion.) Aliasing however 
probably isn't as big of an effect on the BAT PSD's since 
we are averaging over at least 20 points for each time bin. 
Aliasing only becomes a large effect if there are big dif- 
ferences between the exposure of the time bin and the 
gaps between the time bins (i.e. a 1 ksec observation 
every 500 ksec). 

The third complication is the determination of the er- 
rors in the raw periodograms. Papadakis & Lawrence 
(1993) found that only when there are > 20 points av- 
eraged together in a PSD logarithmic bin does the esti- 
mate of the power at a given frequency follow a Gaussian 
distribution. High frequency bins easily meet this crite- 



rion however the low frequency points in the raw PSDs 
contain much less than 20 points so errors cannot be de- 
termined in the usual way by calculating the scatter in a 
bin. 

To account for both windowing effects and the er- 
ror determination, we use a Monte Carlo based tech- 
nique first developed by Uttley, McHardy, & Papadakis 
(2002), hereafter UMP02, called PSRESP that stands for 
Power Spectral Response, and is a method that improved 
upon the response method introduced by Done et al. 
(1992). The basic process is to simulate many light 
curves for each set of model parameters, calculate and 
average together the PSDs of the simulated light curves, 
and determine the goodness of fit for the model parame- 
ters which is defined by a rejection probability. The set 
of parameters with the lowest rejection probability are 
then deemed the best-fit model parameters. In the next 
sections we detail each step of the model fitting. 

3.2.1. Simulating Light Curves 

To simulate the light curves we use the algorithm de- 
scribed in Timmer & Koenig (1995). This algorithm not 
only randomizes the amplitude but also the phase which 
more realistically simulates the noise process that pro- 
duced the observed PSD. The resulting simulated light 
curve is then calculated using the inverse Fourier Trans- 
form. 

The first step is to input a model PSD that has some 



Shimizu, Mushotzky 











Table 2 




















AGN Source Properties 














logLx 


log Lboi 


log Mbh 








logAfn 




Name 


Type 


(ergss-l) 


(ergss i; 


1 (Mo) 


Reference 


Method 


-'jbolZ-^cdd 


(cm-2) 


Reference 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


3C 273 


Blazar 


46.47 


47.62 


8.95±0.09 


P04 


RM 


3.61 


20.5 


B07 


3C 390.3 


Syl 


44.89 


45.86 


8.44±0.14 


WIO 


RM 


0.20 


21.1 


Wi09 


3C 454.3 


Blazar 


47.67 


48.96 


9.18±0.30 


W02 


R-L 


46.08 


20.9 


Wi09 


4U 1344-60 


Sy 1.5 


43.6 


44.42 


7.44±0.10 


VIO 


M-L 


0.07 


22.2 


Wi09 


Cen A 


NLRG 


43.01 


43.76 


7.85±0.24 


C09, N07, K07 


Stars 


0.01 


22.7 


Bll 


Circinus Galaxy 


Sy 2 


42.09 


42.73 


6.23±0.08 


G03 


Maser 


0.02 


24.6 


A12 


Cygnus A 


NLRG 


45.01 


45.99 


9.40±0.12 


T03 


Gas 


0.03 


23.0 


Wi09 


ESQ 506-G027 


Sy 2 


44.12 


45.00 


8.02±0.10 


VIO 


M-L 


0.07 


23.9 


Wi09 


IC 4329A 


Syl 


44.22 


45.11 


8.34±0.21 


099, T02 


M-cr 


0.05 


21.8 


Wi09 


IGR J21247+5058 


BLRG 


44.25 


45.14 


6.58±0.07 


WilO 


R-L 


2.80 


22.1 


Wi09 


MCG-05-23-016 


Sy 1.9 


43.5 


44.30 


6.29±0.70 


WZ07 


Gas 


0.80 


22.2 


Wi09 


MCG+08-11-011 


Sy 1.5 


44.09 


44.96 


8.07±0.02 


WilO 


R-L 


0.06 


21.4 


Wi09 


Mrk 110 


Syl 


44.21 


45.10 


7.38±0.14 


WIO 


RM 


0.40 


20.3 


Wi09 


Mrk3 


Sy 2 


43.74 


44.57 


7.72±0.10 


VIO 


M-L 


0.05 


22.1 


Wi09 


Mrk 348 


Sy 2 


43.91 


44.76 


7.41±0.10 


VIO 


M-L 


0.17 


23.2 


Wi09 


Mrk 421 


Blazar 


44.46 


45.38 


8.56±0.18 


W08 


M-cr 


0.05 


20.9 


Wi09 


Mrk 6 


Sy 1.5 


43.69 


44.52 


8.09±0.02 


WilO 


R-L 


0.02 


22.5 


Wi09 


Mrk 926 


Sy 1.5 


44.76 


45.71 


8.36±0.02 


WilO 


R-L 


0.17 


21.1 


Bll 


NGC 2110 


Sy2 


43.6 


44.42 


7.40±0.10 


VIO 


R-L 


0.08 


22.5 


Wi09 


NGC 3227 


Sy 1.5 


42.57 


43.27 


7.60±0.24 


WIO 


RM 


0.00 


22.2 


Wi09 


NGC 3516 


Sy 1.5 


43.33 


44.12 


7.61±0.18 


WIO 


RM 


0.02 


21.5 


Wi09 


NGC 3783 


Sy 1.5 


43.6 


44.42 


7.45±0.13 


WIO 


RM 


0.07 


21.8 


Wi09 


NGC 4151 


Sy 1.5 


43.11 


43.87 


7.64±0.11 


WIO 


RM 


0.01 


22.7 


Wi09 


NGC 4388 


Sy2 


43.64 


44.46 


7.07±0.10 


VIO 


M-L 


0.19 


23.6 


Wi09 


NGC 4945 


Sy2 


42.37 


43.04 


6.15± 0.16 


G97 


Maser 


0.06 


24.3 


A12 


NGC 5252 


Sy 2 


44.12 


45.00 


8.98±0.21 


C05 


Gas 


0.01 


22.6 


Wi09 


NGC 5548 


Sy 1.5 


43.69 


44.52 


7.89±0.10 


WIO 


RM 


0.03 


20.8 


Wi09 


NGC 6814 


Sy 2 


43.24 


44.01 


7.25±0.12 


WIO 


RM 


0.04 


20.8 


Bll 


NGC 4507 


Sy 2 


43.77 


44.61 


7.70±0.10 


VIO 


M-L 


0.06 


23.5 


Wi09 


NGC 7172 


Sy2 


43.46 


44.26 


7.36±0.10 


VIO 


M-L 


0.06 


22.9 


Wi09 



References. — P04: Peterson et al. (2004); WIO: Woo et al. (2010); W02: Woo & Urry (2002); VIO: Vasudevan et al. (2010); C09: 
Cappollari ot al. (2009); N07: Neumayer et al. (2007); K07: Krajnovic et al. (2007); G03: Greenhill et al. (2003); T03: Tadhunter et al. 
(2003); 099: Oliva ct al. (1999); T02: Tremaine et al. (2002); WZ07: Wang & Zhang (2007); W08: Wagner (2008); G97: Greenhih et al. 
(1997); COS: Capotti ot al. (2005); B07: Bockmann et al. (2007); Wi09: Winter et al. (2009); Bll: Burlon et al. (2011); A12: Ajello et al. 
(2012) 

Note. — Col.(l): Source name. Col. (2): AGN typo taken from HEASARC. Col. (3): 14-195 koV x-ray luminosity from HEASARC. 
Col. (4): Bolonietric luminosity derived using Equation 2. Col. (5): Black hole mass. Col. (6): Reference for black hole mass. Col. (7): 
Method used to determine black hole mass. RM — reverberation mapping, R-L — BLR kinematics using the -Rblr - optical luminosity 
relation, M-L — mass - bulge luminosity relation. Stars — stellar dynamics, Maser — maser dynamics. Gas — gas dynamics, M-cr — mass - 
stellar velocity dispersion relation. Col. (8): Ratio of bolomctric luminosity to Eddington luminosity. Eddington luminosity is determined 
using Equation 3. Col. (9): Column density. Col. (10): Reference for column density. 



continuous shape that wc wish to fit to the observed PSD. 
The normahzation in the model is arbitrary since it wiU 
just multiply through and can be determined later in the 
process. We must also choose for the simulated light 
curve a time resolution, ATgim and overall length, Tgim. 
To fully account for aliasing, we chose ATgim = 800s since 
this is the average exposure time for the raw BAT light 
curves. Any variability from timescales less than 800s 
are averaged over and won't contribute to aliasing. For 
the overall length of the simulated light curve we chose 
lOOOTobs- In this way the simulated long light curve can 
be broken up into 1000 simulated short light curves that 
are Tobs long. Each short light curve will have taken 
into account red-noise leak since it was originally pro- 
duced from a PSD with nonzero power at frequencies 
lower (Tobs)"^- 

Each short simulated curve is then resampled using the 
following procedure. Only those simulated light curve 
points that correspond to a time that was sampled in the 
unbinned observed light curve are used. All others that 
don't are thrown away. This produces a simulated light 
curve that best matches the observed light curve's sam- 



pling so all possible effects including gaps can be included 
modeled. The resampled simulated light curves arc then 
rebinned and gaps are linearly interpolated through in 
the same way as the observed light curve. The PSD for 
each simulated light curve is calculated using equations 1 
&2. 



3.2.2. Determining Goodness of Fit 

After constructing the simulated PSD's, the logarithms 
of them are determined (Psim) and averaged together 
to form Pavg that represents our distorted model PSD 
and has all the effects of windowing and interpolation 
included. The errors in the power at each frequency, 
APsiiii, are equal to the rms spread in the set of Psim- 
Pnoise, which was calculated from Equation 3 and is the 
linear noise power, is then added in linear space as a 



constant to P? 



avg- 



Pavg-l-noise — log (10 



Pn, 



(5) 
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Cen A 
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Figure 3. Top: Distribution of S/N for the sample selected for 
PSD fitting (solid line) and the sample not selected (dashed line). 
Bottom; Distribution of luminosity for the sample selected for PSD 
fitting (solid line) and the sample not selected (dashed line). 
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Figure 4. Observed PSD's of Cen A and NGC 4151 constructed 
using Equation 2 and the BAT lightcurves. Solid line is the PSD 
while the dashed line is the estimated white noise from Equation 3. 
The units on the power are (rms/mean)^(Hz)~-'^ 



The test statistic, Xdisf that is used to compare Pobs 
with Pavg from our model is defined in UMP02 as 



Adist 



E 



(P. 



avg+noise 



(Vi)- Pobs{Vz)) 



AP? 



(6) 



where Vi are the binned frequencies the PSDs are defined 
at. The best-fitting normalization for the model is then 
determined by minimizing Xdist- 

Although Xdist 1^ defined in exactly the same way as the 
traditional x^ statistic, it does not follow the same dis- 
tribution. In fact Xdist 's distribution changes depending 
on the input model parameters. Therefore to determine 
the goodness of fit we use the 1000 Psim to calculate 1000 
Xdist sim which is defined the same as equation 6 except 
that Pobs is replaced by Psim- Each Xdist sim ^^ minimized 
by finding the best-fit normalization in the same way as 
for the observed PSD. The goodness of fit is represented 
by the rejection probability, P, which is the percentile of 
Xdist sim that are less than Xdist ■ P ^^^ ^e thought of as 
analogous to a p-value for a test statistic that measures 
the probability that the test statistic occurred by chance. 
A high P (i.e. low p-valuc) signifies that the model does 
not correctly describe the observed PSD. For this paper 
however we will report Acceptance Probability which is 
just 1 - P. (See UMP02 for a fuU description of the 
fitting procedure.) 



4.1. 



4. RESULTS 
Unbroken Power Law 



For each source we first attempted to fit an unbroken 
power law model which has the form 



PH 



.^0 



(7) 



Aq is the normalization and we chose vq to be equal to 
10^^ Hz to be consistent with previous PSD studies of 
AGN. a is the power law slope which we varied from 
a = 0— 3.0 with a step size of 0.05, because none of the 
observed PSDs show any evidence for a positive slope or 
a slope greater than 3. The best fit parameters for each of 
the AGN are listed in Table 3 and plots of the best fitting 
models with the observed PSDs can be seen in Figures 5. 
Errors on a were determined by calculating the amount a 
needs to change to increase P by la assuming a Gaussian 
distribution. 

CT = V2erf-^(P) (8) 



Rr 



= erf 



1 



V2 



(9) 



The confidence interval can then be determined by 
where P(a) = Pcrit- This method was first used by 
Markowitz et al. (2003b), and it has been shown that 
the errors determined in this way correspond to a 99% 
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Table 3 

PSD Best Fit Parameters for Unbroken Model 



Name 
(1) 


(2) 


log^O 

(Hz-i) 

(3) 


(4) 



3C 273 
3C 390.3 
3C 454.3 
4U 1344-60 
Cen A 

Circinus Galaxy 
Cygnus A 
ESO 506-G027 
IC 4329A 
IGR J21247+5058 
MCG-05-23-016 
MCG-l-08- 11-011 
Mrk 110 
Mrk 3 
Mrk 348 
Mrk 421 
Mrk 6 
Mrk 926 
NGC 2110 
NGC 3227 
NGC 3516 
NGC 3783 
NGC 4151 
NGC 4388 
NGC 4945 
NGC 5252 
NGC 5548 
NGC 6814 
NGC 4507 
NGC 7172 



1 otr + 0.35 
J^-'->3_0.4 
r) 7^ + 0. 65 
'J- '0-0.55 



i.io; 

0.80 
1.35 
0.40 
0.30 
2.40 
0.50 
0.50 
0.15 



0.3 

+0.6 

0.65 

+0.35 

+0.35 

,^0.35 

oo 

0.35 



0.8 

+0.35 

0.35 

0.25 

oo 
+0.4 



1 10+°'' 
0.50+'L'''^ 



0.40: 
0.75 
0.85 
0.701 

2.45^ 



+0.3 
— oo 
+0.3 

0.35 
+0.25 
0.25 
+0.7 
-0.5 
-oo 
-0.85 
0.70+0-2'^ 

0-80lo:35 

0.351°^* 

0.50t°^3 

1.10: 

0.80 

0.95; 

1.10 

0.95 

0.85 

0.70 

0.50 



+0.45 

0.25 
+0.35 

0.3 
+0.25 

0.35 
+0.35 

0.4 

oo 

0.45 
+0.65 
— 00 
+0.3 

0.4 
+0.35 



4.06±0.05 
4.53±0.04 
4.73±0.04 
4.71±0.04 
3.82±0.06 
4.38±0.04 
5.20±0.05 
2.06±0.91 
4.52±0.05 
4.97±0.04 
4.34±0.05 
4.32±0.04 
5.15±0.04 
5.16±0.04 
4.85±0.04 
5.56±0.04 
5.20±0.04 
1.06±1.03 
4.69±0.04 
5.20±0.05 
5.00±0.05 
4.99±0.05 
4.13±0.04 
4.59±0.04 
5.00±0.04 
4.95±0.04 
4.53±0.04 
5.23±0.18 
4.86±0.04 
5.27±0.05 



0.170 
0.728 
0.064 
0.979 
0.379 
0.897 
0.563 
0.164 
0.009 
0.810 
0.936 
0.365 
0.626 
0.956 
0.473 
0.920 
0.986 
0.146 
0.674 
0.292 
0.995 
0.497 
0.205 
0.874 
0.406 
0.544 
0.241 
0.993 
0.912 
0.967 



Note. — Col. (1): Name ot AGN. Col. (2): Best 
fitting power law slope with Ict errors determined using 
the method described in Section 4. Col. (3): Best fitting 
normalization. Col. (4): Acceptance probability for the 
best fit. High values represent a good fit. 



confidence interval when an unbroken power law model is 
used and a 90% confidence interval when a singly broken 
power law model is used (Markowitz 2010, Appendix B). 
The errors on ^g were calculated from the rms spread 
of the 1000 best fitting normalizations found when pro- 
ducing the Xdist distribution. Also instead of showing 
rejection probability, we show the acceptance probabil- 
ity which is just 1 — i?. 

For 3C 454.3 two peaks in the Acceptance Probability 
occured, one at a = 1.1 and the other at 2.75. Both had 
very nearly the same Acceptance Probability (.064 vs. 
.065 respectively) so we chose to use the a = 1.1 fit based 
on visual inspection of the fits. It is well known that 
large values of a can produce large values of Acceptance 
Probability because red-noise leak will cause bigger er- 
ror bars, and thus smaller values of x^jg^. (Vaughan et al. 
2003b; Vaughan & Fabian 2003). The a = 2.75 fit is 
mainly dominated by Pnoiso and the values of the dis- 
torted model don't correspond well with the observed 
PSD values. Figure 6 shows the two peaks in Accep- 
tance Probability and the fits using both peaks. 



We have identified 4 sources which represented "bad" 
fits because most of the power in the fitted model is due 
to Pnoisc shown by a severe flattening in the PSD. These 
are ESO 506-G027, MCG-05-23-016, IC 4329A, and Mrk 
926. For ESO 506-G027, IC 4329A, and Mrk 926 it is 
clear from the figures the fit is bad but for MCG-05- 
23-016 a closer look was needed. Figure 7 shows the 
observed PSD, model fit, and model fit removing any 
distortions, including window effects and Pnoisc- From 
this, we can see that the fit is entirely dominated by 
-Pnoise bccause the undistorted model is below Pnoisc so 
we remove this source from any analysis. The reason 
for these bad fits is either a lack of intrinsic variability 
which is suspected for IC 4329A and MCG-05-23-016, or 
a high Pnoisc which is suspected for ESO 506-G027 and 
Mrk 926. Pnoisc may not be the same in the BAT data 
due to the method via which the catalog is constructed 
(Tueller et al. 2010) which involves cleaning the data. 

24 out of 26 sources (our original 30 sources minus the 4 
with bad fits) had Acceptance Probabilities greater than 
20% indicating the unbroken model is a sufficient model 
to characterize the hard X-ray PSDs of AGN. The two 
sources with less than 20% acceptance are 3C 273 and 
3C 454.3 which will later be tested with a singly broken 
model. 

The average power law slope for all AGN with a best 
fit slope was 0.78±0.29 where the error is the standard 
deviation of the slopes. If we restrict the average to only 
those slopes that were fully constrained then it changes to 
0.92±0.22. However this might be biased against flatter 
slopes because those are more likely to be unconstrained. 

16 out of 26 sources had their PSD power law slopes 
fully constrained. 10 of the other 14 sources only re- 
sulted in upper limits on the slope, and the other 4 un- 
constrained sources only resulted in lower limits on the 
slope. These are indicated in Table 3 by a ±cx) in one of 
the error bars. Flat slopes are difficult to constrain on 
the lower end because the effect of Pnoisc is to flatten the 
PSD out. So if the observed PSD power is near the level 
of Pnoise, many different slopes will provide a good fit. 

4.2. Singly Broken Power Law 

For 4 AGN (3C 273, Cen A, NGC 4151, and 3C 454.3) 
wc tested a singly broken power law model with the form 






A, 



(10) 



ly > vt,. 



where Ai is the normalization, vi, the frequency at which 
the PSD changes from a high frequency power law slope 
of — /3 to a low frequency slope of —a. We chose these 
AGN due to their steep best fit slope from the unbro- 
ken model and/or their low Acceptance Probability. We 
suspected that possibly our low frequency resolution for 
the PSD could be smoothing out a potential break fre- 
quency thus giving an "average" best fit slope somewhere 
between 1 and 2 depending on where the break occurs. 

Figures 8, 9, 10, 11, and Table 4 show the results 
from this model. Fitting the unbroken power law model 
took on average 1 day per source on a machine with 3.0 
GHz processing speed and 8 cores. Using a singly bro- 
ken power law model with 2 more free parameters could 
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Figure 5. Observed PSD's (solid line) plotted with best fitting distorted model PSD (dots with error bars). The distorted model contains 
all windowing effects and error bars are Icr uncertainties. Dashed line indicates estimated white noise level for each AGN due to measurement 
error. See Table 3 for a list of best fitting parameters and acceptance probabilities. 
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Figure 5. continued 
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Figure 5. continued 
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Figure 5. continued 
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Table 4 

PSD Best Fit Parameters for Singly Broken Model 



Name 
(1) 


Type 

(2) 


a 
(3) 


logi^b 
(Hz) 

(4) 


logAi 

(Hz-i) 

(5) 


\-R 

(6) 


Cen A 


NLRG 


i-3ir3 


a n+oo 

-6-0_o.4 


3.97±0.05 


0.262 


3C 273 


Blazar 


0.6t«J 


-6.8t°i 


5.55±0.02 


0.680 


NGC 4151 


Sy 


1.5 


1 1+0.5 
J^-^-0.3 


-6.0tS°3 


4.26±0.04 


0.203 


3C 454.3 


Blazar 


1 2+°° 


-6.4±CS 


5.27±0.05 


0.096 


3C 390.3 


Blazar 


1.0 


-6.2t-2 


4.64±0.04 


0.621 


4U 1344-60 


Sy 


1.5 


1.0 


-6.otr=4 


4.64±0.04 


0.802 


Circinus Galaxy 


Sy 


2 


1.0 


-6.2l- 


4.64±0.04 


0.621 


Cygnus A 


NLRG 


1.0 


-8.0+^ 


7.70±0.13 


0.000 


ESQ 506-G027 


Sy 


2 


1.0 


-8.0l^ 


7.46±0.13 


0.030 


IC 4329A 


Sy 


1 


1.0 


-6 2+°-^ 

"•^-0.1 


4.84±0.04 


0.001 


IGR J21247+5058 


BLRG 


1.0 


-8.0t- 


7.58±0.13 


0.000 


MCG-05-23-016 


Sy 


1.9 


1.0 


-8.0t- 


6.40±0.13 


0.000 


MCG+08-11-011 


Sy 


1.9 


1.0 


-6-5-0.9 


5.10±0.05 


0.453 


Mrk 110 


Sy 


1 


1.0 


-6-0+^ 


4.91±0.04 


0.073 


Mrk3 


Sy 


2 


1.0 


-8.0l^ 


7.57±0.13 


0.000 


Mrk 348 


Sy 


2 


1.0 


-6.0tS°5 


4.88±0.04 


0.073 


Mrk 421 


Blazar 


1.0 


-6.0tS°3 


5.66±0.05 


0.337 


Mrk 6 


Sy 


1.5 


1.0 


-e.otrg 


5.06±0.04 


0.786 


Mrk 926 


Sy 


1.5 


1.0 


-6.4l^ 


4.83±0.04 


0.027 


NGC 2110 


Sy 


2 


1.0 


-6.0tS°4 


4.76±0.04 


0.085 


NGC 3227 


Sy 


1.5 


1.0 


-6.0tS°5 


5.22±0.04 


0.124 


NGC 3516 


Sy 


1.5 


1.0 


-e.otai 


5.05±0.04 


0.001 


NGC 3783 


Sy 


1.5 


1.0 


-6 0+°° 


5.04±0.04 


0.003 


NGC 4388 
NGC 4507 


Sy 
Sy 


2 
2 


1.0 
1.0 


-6.0ta4 
-6-0tS°5 


4.65±0.04 
4.88±0.04 


0.444 
0.080 


NGC 4945 
NGC 5252 


Sy 
Sy 


2 
2 


1.0 
1.0 


-6 0+°° 
0-'J_o.4 

-6.4ta5 


5.13±0.04 
5.60±0.05 


0.263 
0.769 


NGC 5548 


Sy 


1.5 


1.0 


-6.5t^ 


5.12±0.05 


0.258 


NGC 6814 


Sy 


2 


1.0 


-6-il:?°5 


5.38±0.18 


0.954 


NGC 7172 


Sy 


2 


1.0 


6 1+°° 


5.46±0.04 


0.005 



Note. — Results from fitting a singly broken model. For Cen A, NGC 
4151, 3C 273, and 3C 454.3. ct, the power law slope below the break fre- 
quency and i^b, the break frequency were freely allowed to vary while the 
high frequency slope was fixed at 2. For the rest of the sources a was fixed 
at 1.0 and only i^^ was allowed to vary. 
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Figure 6. Fit results for 3C 454.3. Top: Curve of Aeceptance 
Probability against power law slope (a). Notice the two peaks in 
the curve at 1.1 and 2.75. Bottom: Comparison of the fit with 
a = 1.1 (squares, colored red in online version) and a = 2.75 
(circles, colored blue in online version). The black solid curve is 
the observed PSD and the dashed line is the estimated P„oise level. 




-7.0 -6.5 

log (y) 



Figure 8. Fit results for 3C 273 using the singly broken model 
described by Equation 10. Top: Contours of constant Acceptance 
Probability for the different pairs of break frequency (uf,) and low 
frequency slope (a). 3 contours were drawn at 32, 10, and 1% 
the peak of Acceptance Probability. Stars indicate the location of 
the best fit parameters. Bottom: Observed PSD (solid line) plot- 
ted along with best fit distorted model (dots with error bars) and 
undistorted model (dot-dash). Horizontal dashed line represents 
the Poisson noise level. 
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take on the order of months, so to reduce the computa- 
tional time, we froze /3 ~ 2 reducing the number of free 
parameters to 2, vt and a. Freezing the high frequency 
power law slope seems to be a valid decision since this 
is the typical value measured for the PSD in the 2-10 
keV band (Gonzalez-Martin & Vaughan 2012) ,and we 



are most likely to be sensitive to the low frequency slope 
rather than the high frequency one since we are prob- 
ing the longest timescales. Vf, was stepped through ev- 
ery 0.1 in the log between -8 and -6, approximately the 
lowest and highest temporal frequencies sampled in our 
observed PSDs, and a was stepped through every 0.1 (as 
opposed to 0.05 with the unbroken model) between and 
3. 

Both Cen A and NGC 4151 resulted in essentially the 
same fit as for the unbroken model. The best fit a's 
were exactly the same as for the unbroken model (since 
we stepped through a with stepsize 0.1, 1.35 was not 
an option for Cen A) as well as the normalization. The 
best fit break frequencies also occured at the edge of the 
parameter space demonstrating that the break frequency 
must occur at a frequency outside the range sampled by 
the PSD's. 

The singly broken model only gave marginal improve- 
ment for 3C 454.3 as well (9.6% vs. 6.4% Acceptance 
Probability). While a best fit break in the PSD was mea- 
sured at log(i^b) = —6.4, Fig 11 (left) clearly shows that 
this is simply where the Poisson noise begins to domi- 
nate as evidenced by the intersection of the dot-dashed 
line and the dashed line. Also the break was completely 
unconstrained with the 90% confidence interval spanning 
the entire range of the parameter space. We also looked 
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Figure 9. Same as Figure 8 but for Cen A 




Figure 10. Same as Figure 8 but for NGC 4151 



at the second peak in the Acceptance ProbabiUty which 
occurred at \og{i'b) = ~7.4 and a = 0.2 shown in the 
right column of Fig 11. The Acceptance Probabihty for 
this peak was 7.2% so still an improvement from the un- 
broken model. The low frequency slope found for this fit 
would be much flatter than previously seen and could be 
analogous to the low frequency breaks seen in Galactic 
Black Holes (GBH) while in their low-hard state. How- 
ever this break is only based on the 2 or 3 lowest fre- 
quency PSD points which are also the points that are 
the least statistically significant since they represent the 
longest timescales. Due to this and the low Acceptance 
Probability we conclude that there is no break in 3C 
454.3. 

3C 273 on the other hand resulted in a much bet- 
ter fit with the singly broken model with an acceptance 
probability of 68% compared to 17% for the unbroken 
model. The best fit break frequency and low frequency 
slope occurred at log(i^;,) = —6.8 and a = 0.6 respec- 
tively. This corresponds to a timescale of 73 days or 
about 2.5 months. The break was well constrained be- 
tween log(i'h) = —6.5 and -7.2 but the slope was only 
constrained to an upper limit of 1.3 which is still flatter 
than the slope found for the unbroken model at 1.35. 

For the rest of the sample, we also fit a singly broken 
model but only allowed log(i/(,) to vary between -8.0 and 
-6.0 at 0.1 steps to determine an approximate lower limit 
for the break frequency and fixed a = 1.0. The results of 
these fits arc also displayed in Table 4. 7 of the sources 
resulted in completely unconstrained break frequencies. 
Most of these sources were either the 4 sources we flagged 
as bad fits for the unbroken model, or had significantly 
flat slopes {a < 0.5) so neither 1.0 nor 2.0 was a close 
approximate to the slope of the PSD. 

We estimate the lower limit of the break frequency on 
our sample to be log(i/t) = —6.35. This was found by av- 
eraging the upper limits of the break frequencies for all of 
the sources where a lower limit was constrained including 
Cen A and NGC 4151. This corresponds to a timescale of 
26 days and is very close to the upper edge of our sampled 
frequency band where Pnoise usually dominates. 

4.3. Energy Dependence of PSD 

For 5 sources (Cen A, NGC 4151, Mrk 421, NGC 3227, 
and NGC 2110), the S/N was high enough to split into 
three separate energy bands and determine best fit pa- 
rameters for the PSD of each energy band. Table 5 lists 
the best fitting power law slope, normalization and Ac- 
ceptance Probability for each source in each band. We 
split the 14-150 keV fufl band into 14-24 keV (Band 1), 
24-50 keV (Band 2), and 50-150 keV (Band 3) bands. 
These represented the best bands to get nearly equal 
S/N in each band but also allows for possible separation 
of spectral components because refiection should approx- 
imately be most important in the 24-50 keV band. Fig- 
ure 12 plots the best fit models for each source in each 
band. 

All of the AGN were fit very well with an unbroken 
power law again, except for NGC 3227 in the 50-150 keV 
band but this was due to poor S/N in the PSD rather 
than an ill-fitting model. The average power law slopes 
for each band are 0.88±0.24, 1.12±0.29, and 0.94±0.33 
for Band 1,2, and 3 respectively with errors the standard 
deviation of the sample. There seems to be some evi- 
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Table 5 

Best Fit Parameters: Energy Dependence 



14 - 24 keV 



24 - 50 kcV 



50 - 150 keV 



Name 
(1) 



log^O 
a (Hz-i) 1-R 

(2) (3) (4) 



logAo 
a (Hz-i) 1-R 

(5) (6) (7) 



log^O 
a (Hz-i) 1~R 

(8) (9) (10) 



Cen A 
NGC 4151 
Mrk 421 
NGC 3227 
NGC 2110 



1.15 
1.10 
0.95 
0.65 
0.55 



+0.35 

0.5 
+0.25 

0.3 
+0.25 

0.25 
+0.35 

0.4 
+0.3 

0.3 



4.17±0.04 

4.35±0.04 
5.43±0.04 
5.47±0.05 
5.00±0.05 



0.317 
0.577 
0.943 
0.735 
0.675 



1.55 
1.30 
0.80 
0.80 
1.15 



+0.35 
0.25 
-0.3 
-0.3 

+0.3 
0.25 

+0.55 
0.3 

+0.4 
0.45 



3.54±0.08 
3.99±0.05 
5.64±0.04 
5.18±0.04 
4.46±0.04 



0.399 
0.365 
0.768 
0.306 
0.139 



1.55 
1.05 
0.70 
0.70 
0.70 



+0.65 

-0.35 

+0.4 

0.3 

+0.4 

-0.35 

+0.5 

0.3 

+0.35 

0.2 



3.66±0.08 
4.27±0.05 
5.90±0.04 
6.00±0.04 
4.69±0.04 



0.561 
0.817 
0.977 
0.006 
0.674 



Note. — Best fitting parameters for the PSDs in each energy band. 



dence for the increase of the slope between bands 1 and 
2, however it is only significant at the la level. Moreover, 
all three bands are consistent with a power law slope of 
1.0 so we conclude there is no detectable dependence of 
the PSD with energy. 

4.4. Correlation with BH Mass, X-ray Luminosity, 
Accretion Rate, and Column Density 

To determine any significant correlations with prop- 
erties of the AGN, we used the integrated PSD as the 
dependent variable rather than any of the parameters of 
the PSD itself. Integrating over the best-fit PSD gives 
us an estimate of the excess variance produced by the 
source itself, in other words variance above that expected 
by measurement error and noise. Excess variance is then 
defined as 



^2 



P{v)dv 



(11) 



which can be analytically integrated with an unbroken 
power law model of the form in Equation 7 to find 



2 



An 



(1 



(l-a) 



(12) 



For all of the AGN we substituted in for A^ and a 
the best fit normalization and power law slope found in 
Table 3 however we excluded from the analysis the 4 
AGN discussed previously that had bad fits to the PSD. 
To test that we were getting reasonable values of the 
excess variance, we determined it using the time domain 
definition (Nandra et al. 1997; Edelson et al. 2002). 
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Figure 12. Results from fitting tlie 14-24 kcV, 24-50 keV, and 50-150 kcV PSDs for Cen A, NGC 4151, Mrk 421, NGC 3227, NGC 
2110. Circles and solid lines correspond to the 14-24 band (colored red in the online version), squares and dashed lines correspond to the 
24—50 keV band (colored blue in the online version), and the triangles and dot-dashed lines correspond to the 50-150 keV band. The non 
horizontal lines plot the observed PSD, the straight horizontal ones the Poisson noise level, and the markers with error bars are the best 
fit model PSD. 



2 



52- 



^I- 



where S'^ is the sample variance 



TV 



S^ 



N- 



tE( 



M)^ 



(13) 



(14) 



and cr^j.j. and fi are defined as before. The errors on 
the integrated PSD variance were determined by prop- 
agating the errors on the normahzation and the slope 
through Equation 12 while also taking into account the 
anti-correlation between normalization and power law 
slope. The errors on the time domain variance were cal- 



culated using the following equation 



err{a- 



xs) 



\ 





2F 



(15) 

where N is the number of light curve points, and i^var 
is the square root of a^g (Vaughan et al. 2003a). But 
this only takes into account the error due to measure- 
ment uncertainties. There is also uncertainty due to the 
random scatter intrinsic to the red noise process that 
produces the variability. For this uncertainty we used 
Table 1 from Vaughan et al. (2003a). Our light curves 
contained 400 points and the PSDs had a slope of about 
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Table 6 

Correlation between Hard 

X-Ray Variability and Source 

Properties 



Variable 


Pa 


Po 


(1) 


(2) 


(3) 


All Sample 






LboI 


0.14 


0.49 


Mbh 


0.08 


0.70 


LBol/LEdd 


-0.10 


0.63 


Nh 


-0.16 


0.44 


Seyferts Only 






LboI 


-0.01 


0.97 


Mbh 


-0.04 


0.87 


LBol/LEdd 


-0.23 


0.32 


Nh 


-0.09 


0.71 



Note. — Col. (1): Source 
property used in correlation anal- 
ysis with hard X-ray variance cal- 
culated by integrating the best fit 
PSD. Col. (2): ps is the Spear- 
man rank correlation coefficient 
Col. (3): Pp the null probability. 
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Figure 13. Ratio of the excess variance determined by integrating 
the best-fit PSD and by using the variance of the light curves. 
Dashed line indicates a ratio of 1. 



1, SO we used an uncertainty of 0.18 which is halfway be- 
tween the uncertainty of a light curve with 100 points and 
one with 1000 points (Vaughan et al. 2003a). We then 
summed in quadrature the uncertainty from Equation 15 
and 0.18 to determine the total uncertainty. Figure 13 
shows the relationship between the PSD integrated and 
time domain cr^g. While there is some scatter in the re- 
lationship, it seems that both indicators agree fairly well 
with each other within the error. 

Table 6 shows the results of a Spearman rank corre- 
lation test between the hard X-ray variance and bolo- 
metric luminosity, black hole mass, accretion rate with 
L^ioi/Lcdd as the indicator, and column density. We 
looked for correlations using both the whole sample of 
AGN and a sample of just Seyfert galaxies because radio- 
loud AGN could have variability introduced due to other 
processes associated with the jet especially in the energies 
studied here. Figures 13 & 14 display the 4 correlations 
for both samples. 

As can be seen in the Figures and Table, there seems 
to be no significant correlation between any of the source 
properties and hard X-ray variance. None of the corre- 



lations have a null hypothesis probability, Pp < 0.3 and 
the largest correlation coefficient is merely 0.23. This 
result agrees with those of Markowitz & Edclson (2004) 
and Matcos et al. (2007) that found the anticorrelation 
between variability amplitude and luminosity and black 
hole mass weakens as longer timescales are sampled. This 
can be explained as a consequence of the simple scaling 
of PSD break timescales with luminosity and black hole 
mass. On short timescales, the difference in PSD breaks 
between different AGN will manifest itself when deter- 
mining variability amplitude since this is merely an inte- 
gration of the PSD. However as one moves to longer and 
longer timescales away from the PSD break, variability 
amplitude increases (as it should for red noise processes) 
but eventually saturates, and the dependence on lumi- 
nosity and black hole mass disappears. This makes sense 
in our sample as well since we report only one object 
where the PSD break is in the frequency band sampled 
by the BAT light curves. 

5. DISCUSSION 
5.1. Comparison with 2-10 keV Variability 

We cross-referenced our sample with the samples 
of previous PSD studies in the 2-10 keV band and 
found 9 AGN's PSD had previously been measured. 
Table 7 lists these AGN and their PSD parame- 
ters. Most of the 2-10 keV studies used RXTE light 
curves to determine low frequency power and either 
XMM or Chandra light curves for the high frequency 
power. Gonzalez-Martin & Vaughan (2012) (GM12) 
used shorter length XMM light curves so the observed 
PSD only covered the high frequencies, and when the 
authors fit the PSD, they fixed the low frequency slope 
to a value of -1.0. Figure 16 shows the PSDs measured in 
the two bands for the timescales sampled in this study. 

In general, there does not seem to be a significant dif- 
ference in the power law slopes of the soft and hard 
X-ray PSDs. We find an average slope of 0.8, sim- 
ilar to the low frequency slope of '--^l found in the 
lower energy studies (Uttley et al. 2002; Markowitz et al. 
2003b; McHardy et al. 2004; Uttley & McHardy 2005; 
Markowitz et al. 2007; Markowitz 2009). Kelly et al. 
(2011), using Bayesian analysis and modeling light curves 
as a linear combination of Gaussian Ornstein-Uhlenbeck 
processes determined the slopes and break frequencies of 
optical and soft X-ray light curves of AGN. This mod- 
eling technique allows for a better interpretation of the 
actual physical processes occurring in the AGN. They 
also find a power law slope of 1.0 describes very well the 
long timescale PSD in the 2-10 keV energy range. This 
suggests the long timescale variability process is the same 
between the 2-10 and 14-150 keV energy bands. 

We also compared the intrinsic variances of the AGN 
in both energy bands, again using the excess variance 
calculated directly from the light curves to get a better 
estimate of the uncertainty. For the 2-10 keV PSDs, we 
integrated the PSDs we found in the literature, taking 
into account the break frequency when needed. Some of 
the AGN were best fit with a slowly bending power law 
model instead of a sharp break, however when integrat- 
ing the PSDs, we strictly used the sharp break model to 
get a fairly close estimate, especially since the break is 
usually near the edge of our frequency range or outside 
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Figure 14. Correlation between excess variance of the whole sample BAT AGN and their Nh (top left), Lgoi (top right), Mbh (bottom 
left), and LboiI LEdd (bottom right). Circles (colored blue in online version) are Seyfert 1-1.5, diamonds (colored red in online version) are 
Seyfort 1.8-2, stars (colored black in online version) are Blazars, triangles (colored cyan in online version) are BLRG, and squares (colored 
magenta in online version) are NLRG. See Table 6 for results of statistical analysis. 
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it completely. 

Figure 17 shows a comparison between the variances of 
the two energy bands with a dashed line indicating equal 
variances. All but two sources (NGC 4151 and Cen A) 
show an apparent increase in variance when the energy 



band increases from 2-10 keV to 14-150 keV. This result 
agrees with that found by Caballero-Garcia et al. (2012) 
which only looked at the brightest AGN in the BAT sam- 
ple and is surprising because variance had been seen pre- 
viously to decrease as energy increases (McHardy et al. 
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Table 7 

2-10 kcV PSD Parameters 









2-10 keV 








14-150 keV 




Name 


a2-io 


iogi^i, 


/^2-ia 


Ai 


Hetcrencc 


"14-150 


Ao or Ai 


lOgt/fc 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


3C 273 


1.0 


-6.2 


2.4 


10* 


M04 


0.6l°J 


3.16 X 10^ 


-6.8 4 


3C 390.3 


1.0 


-6 6+°* 

"•"-0.2 


9 4+0.3 


2.51 X 10* 


G09 


0.75t°otl 


3.42 X 10* 






Cen A 


0.91^2 


-6 2+0-3 


9 c:+0.4 
^■"-0.1 


7.90 X 10-* 


Rll 


1 oc+0.35 
J^-'l"_0.35 


6.61 X 10^ 






Circinus Galaxy 


1.0 


o on+0.13 
o.oa_Q j3 


o 90+0.12 


8.10 


GM12 


0.4«;f 


2.39 X 10-* 






NGC 3227 


1 n+0'2 

^^•"-0.4 


_4 Kq+0.36 


> 1.9 


371 


U05 


0-80l°J5 


1.59 X 10^ 






NOG 3516 


^•^-0.3 


c 7+0.4 
"■'-0.3 


9 n+0.55 
^■'^-0.2 


7.90 X 10^ 


Ma03 


0.35l°:|5 


9.96 X 10-* 






NGC 3783 


o-8ol°:| 


t: 9 + 0.88 


9 (5+0.6 


1.10 X 10^ 


S07 


o.5tg;i 


9.83 X 10* 






NGC 4151 


I ;^ + 1.35 


-5 9+°-^ 

"■^-0.7 


9 1+1.9 


1.30 X 10* 


Ma03 


1 1+0.45 

J^-^-0.25 


1.34 X 10* 






NGC 5548 


1 ic;+0-5 


f; 99 + 0.60 
-O.^^_0.40 


2.051H 


2.50 X 10-* 


Ma03 


0.951^:°^ 


1.01x10^ 







References. — M04: Marscher ot al. (2004); G09: Gliozzi ct al. (2009); Rll: Rothschild ot al. (2011); GM12: 
Gonzalez-Martin & Vaughan (2012); U05: Uttloy & McHardy (2005); Ma03: Markowitz et al. (2003a); S07: Summons et al. 
(2007) 

Note. — Comparison between the 2—10 keV PSDs found from a literature search and 14—150 keV PSDs from this work. Col. 
(1): Name of source. Col. (2): Low frequency power law slope for 2-10 keV PSD. Col. (3): Break frequency for 2-10 keV PSD. 
in Hz. Col. (4): High frequency power law slope for 2-10 keV PSD. Col. (5): Normalization for 2-10 keV PSD in Hz"^ Note 
that the normalization here is defined as the variability power at the break frequency, not at 10~ Hz like in this work. Col. (6); 
Reference for the 2-10 keV PSD. Col. (7) and (8): 14-150 keV PSD parameters found in this work. 



2004; Inoue et al. 2011). Using Suzaku light curves, 
Miyazawa et al. (2009) studied the temporal and spectral 
variations of 36 AGN in 3 different energy bands (0.5-2 
keV; 2-10 keV, and 15-50 kcV). Their sample overlaps 
with our sample, and their highest energy band overlaps 
the BAT band. They as weU found that the 15-50 kcV 
variability is much lower than that in the lower energies. 
However all of these studies used much shorter timescales 
than the one presented here. Whereas our timescales 
are on the order of years, previous studies are only on 
scales of less than half a day. This seems to indicate that 
the variability process is timescale and energy dependent, 
where at shorter timescales lower energies display higher 
variability and on longer timescales variability is at least 
equal across all energies. 

For 3C 273 we were able to compare the break frequen- 
cies between the 2-10 keV and 14-150 keV band. In- 
terestingly the break frequency seems to have decreased 
from logi^f, ~ -6.2 to -6.8 (18 days to 73 days). How- 
ever as shown in Figure 17 the total variance possibly in- 
creased or at the least stayed the same. This is probably 
due to the decrease in both the low and high frequency 
slope where the low frequency slope decreased from 1.0 
to 0.6 (although only an upper limit of 1.3 could be con- 
strained) and the high frequency slope decreased from 
2.4 to 2.0 (although this is because we fixed the slope at 
2.0). We caution the reader that without error bars from 
the 2-10 keV study, a valid comparison cannot be made. 

5.2. Model for Variability? 

Lyubarskii (1997) suggests a model for the variabil- 
ity of accreting objects in which inwardly propagating 
fluctuations of the accretion rate couple at each radii 
of the accretion disk and modulate the X-ray emission. 
This model successfully explains the linear rms-flux re- 
lation seen in all GBH's and AGN (Uttlcy & McHardy 
2001) and can explain the energy dependent time lags 
(Kotov et al. 2001). Arevalo & Uttley (2006) further in- 
vestigate this model and determine that if the emissivity 
of the X-ray emitting corona is radially dependent, then 



there should be a change in the shape of the observed 
PSD with energy. The break frequency correlates with 
the radius at which the emissivity drops off, thus for ex- 
ample, if the harder X-rays originate from a more com- 
pact region than softer X-rays, then the break frequency 
should increase as energy increases since there will be 
more variability on shorter timescales due to the small 
radius. 

Our results do not rule out the possibility of this 
model but only due to the fact that we detected only 
one break frequency in 3C 273. The BAT PSDs were 
limited to the longest of timescales with the short- 
est timescale being ~week. This is about the longest 
timescale that a break has been measured in softer 
bands with most breaks occurring on shorter timescales 
(Gonzalez-Martin & Vaughan 2012). The fact that the 
power law slope stays roughly constant around 1.0 does 
agree with the results of Arevalo & Uttley (2006) who 
found a constant low frequency slope with energy. 

3C 273's measured break frequency, while seeming to 
decrease with energy, is still formally consistent between 
the 2-10 and 14-150 kcV band especially given the lack 
of error bars from Marscher et al. (2004). Because of the 
large uncertainty though, the possibility remains that 3C 
273 exhibits different variability mechanisms in the two 
bands, maybe due to its radio jets. 

Ishibashi & Courvoisier (2009) recently proposed an- 
other model for both accretion onto a black hole and 
variability. Instead of a smooth accretion disk flowing 
into the black hole, they use clumps of material falling 
in. These clumps collide and first form optically thick 
shocks that produce the optical/UV photons typically 
seen in AGN and GBH's. The optically thick shocks 
then become optically thin while the whole time elec- 
trons are being heated. The optical/UV photons origi- 
nally emitted in the optically thick shock provide the seed 
photons for inverse Gomptonization and are upscattered 
into X-rays much the same way in the usual accretion 
scenario. Variability however is attributed to the actual 
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Figure 16. Comparison between the PSDs measured in the 2-10 keV band (dashed) and the 14-150 keV band (solid). 



Comptonization process, and the PSD break timescale is 
associated with the Compton coohiig timescale of the 
heated electrons. This model has been shown to re- 
produce the empirically found correlations between PSD 
timescale, black hole mass, and accretion rate through 
simple physics (Ishibashi & Courvoisier 2012). 

However our results with their limited frequency range 
are not able to discern between the two different models. 
What is needed are high quality light curves at these 
high energies on shorter timescales to be able to measure 
break frequencies in a large sample of AGN. 

5.3. Variance vs. Type 

Figure 18 plots the histograms of excess variances 
based on the type of source. As before we grouped all 
Seyfert 1, 1.2, and 1.5 as Seyfert 1 and all Seyfert 1.8, 
1.9, and 2 as Seyfert 2. But we also grouped all Blazars, 
NLRG, and BLRG as radio loud objects. From the his- 



tograms there does not appear to be any difference in 
the variances of the different types of AGN. The one ob- 
ject with the largest variance is Mrk 421 which has been 
known to exhibit extremely strong variability in the form 
of flares over many timescales (Gui 2004). However it is 
interesting to note that even though Mrk 421 is the only 
source where radiation is primarily beamed, the slope of 
the PSD for Mrk 421 (0.8) is not different than the slopes 
of any of the other AGN. it is just the normalization that 
is significantly higher. This suggests that long timescale 
variability might be caused by the same process no mat- 
ter the type of AGN which agrees with the results seen 
in § 4.4. 

5.4. Comparison with GBH 

GBH's due to their closeness and large fluxes have 
much higher quality spectra. They exhibit state tran- 
sitions characterized by their intensity (high or low) and 
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Figure 18. Distribution of the excess variance by type of source: 
Scyfert 1, Scyfcrt 2, and radio loud AGN. We grouped all Seyfert 
1, 1.2, and 1.5 as Seyfert Is and Seyfert 1.8, 1.9, and 2 and Seyfert 
2s. Radio loud objects include Blazars, BLRG, NLRG. 



shape (soft or hard) of the energy spectrum. The differ- 
ent states also show different variability that change the 
PSD (Cui et al. 1997; Pottschmidt et al. 2003). In the 
high/soft state, the PSD is best fitted by a singly-broken 
or bending power law while in the low/hard state it is 
best fitted by a doubly-broken power law. For nearly all 
AGN, in the 2-10 keV band, their PSD are similar to the 
high/soft state of GBH because they are best described 
by a singly broken power law (see McHardy et al. (2004) 
and Markowitz et al. (2003b)). Currently the only AGN 
found to be best fit by a doubly-broken power law 
is Ark 564 (Pounds et al. 2001; Papadakis et al. 2002; 
McHardy et al. 2007). 

Since our PSD's don't show any evidence for a low 
frequency break, they are consistent with the PSD's 
of GBHs in the high/soft state if the PSD is constant 
with energy. This is in agreement with studies in the 
2-10 keV band including Markowitz et al. (2003b) and 
Uttley & McHardy (2005). Previous work looking at 
the energy dependence of Cygnus X-l's high state PSD 
also observed a constant low frequency power law slope 
with a possible increasing break frequency with energy 
(Nowak et al. 1999; McHardy et al. 2004). Our finding 
that the PSD slope is relatively constant with energy 



agrees with this, and although we did not detect any 
break frequencies, the lack of detection indicates a lower 
limit for the break frequency at 10~^. With the excep- 
tion of 3C 273, our study seems to indicate that all of 
our AGN are in the high state since we don't find any 
evidence for a low frequency break that distinguishes be- 
tween the low and high state. Zhang (2011) used archival 
RXTE All Sky Monitor light curves to determine the 
variance at roughly the same timescales as our BAT light 
curves. They also find that the measured variability 
seems to indicate that AGN are scaled up versions of 
only the high states of GBHs. They determine that if 
AGN were in the low state with a low frequency break 
timescale, there would still be an observed correlation be- 
tween the variance and black hole mass and luminosity 
whereas the high state with a singly broken PSD would 
not produce any correlation at long timescales. This is 
exactly what we observe in our sample at higher energies 
using the more rigorous PSD analysis. 

6. CONCLUSIONS 

We have for the first time calculated and fit the PSDs of 
AGN in the hard X-ray (14-150 keV) regime. The PSDs 
were fit very well with an unbroken power law model 
with an average slope of ~^ — 1 and only one object, 3C 
273, was found to require a break in the PSD. The lack 
of break frequencies at the longest of timescales points 
toward a majority of AGN being in the high state seen 
in GBHs. The slopes found are similar to those found 
in 2-10 keV PSD studies and the total variance at these 
long timescales does not show any evidence of changing 
with energy. 

We did not find any correlation between the hard X-ray 
variance and properties of AGN previously seen in lower 
energy variability studies. It seems that whatever pro- 
cess is producing long timescale variability is completely 
independent of the standard measures of AGN including 
luminosity, black hole mass, accretion rate, and column 
density. This also holds when looking at the distribution 
of the strength of variability with AGN classification. 

Unfortunately our hard X-ray variability study has not 
provided much more insight into verifying the processes 
that are producing the variability since only one PSD 
break was detected. However we have provided an upper 
limit of ~ 26 days for the break timescale of the hard 
X-ray PSDs leaving open the possibility that the break 
timescale still decreases with energy (break frequency in- 
creases) . 

Further work needs to be done in extending the fre- 
quency range of the PSDs. With the recent launch of 
NuStar this will be possible because it will be much more 
sensitive than BAT. Using WebPimms we predict for a 
source with flux 2x 10~^^ ergs/s/cm^ a count rate of .07 
cts/s compared to 3x 10~^ cts/s for BAT. This will allow 
us to sample the high temporal frequencies just as XMM 
provided for the 2-10 keV band. Combining the BAT 
and NuStar PSDs would give us a broadband PSD span- 
ning at least 6 decades of frequency allowing for better 
constraints on models for AGN. It will also allow us to 
significantly constrain the break frequency if one exists 
and determine if the same relation between black hole 
mass, accretion rate, and break frequency exists in the 
hard X-ray band as does in the soft X-ray band found by 
McHardy et al. (2006). 
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